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Abstract 

It is found that combining an excitation-energy sum rule with Fetter's trial 
wave function gives almost exact low-lying collective-mode frequencies of a 
trapped Bose-Einstein condensate at zero temperature. 

Realizations of Bose-Einstein condensation (BEC) in trapped atomic gase have en- 
abled us to study many-body physics of weakly interacting bosons with unprecedented pre- 
cision. Experiments on elementary excitations of BEdi present a new challenge because 
inhomogeneity and rotational symmetries of the system must be explicitly taken into ac- 
count. A number of theoretical studies have been reported in literature JUll The results 
of numerical analyses^! based on the Bogoliubov approximation@'El have so far found the 
best agreement with those of experiments-^ Stringari^l obtained exact analytic expressions 
of collective-mode frequencies in the Thomas-Fermi (TF) limit.&B Perez-Garcfa et al.0 as- 
sumed a Gaussian trial wave function which yields correct collective-mode frequencies at the 
weak and strong coupling limits but overestimates them at an intermediate regime. In this 
Letter, it is found that combining an excitation-energy sum rule!! with Fetter's trial wave 
function^ gives the almost exact frequencies of a trapped Bose system at zero temperature. 
Our method will be shown to yield collective-mode frequencies in excellent agreement with 
experimental results as well as numerical ones based on the Bogoliubov approximation. 
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Consider a Bose system described by the Hamiltonian, H = T + U + H int , where T = 
-(h 2 /2M) E l V? , ?7 = (M/2) Ei(^ + u 2 y 2 + cu 2 z zf) and H int = (27rh 2 a/M) - r,-) 

describe the kinetic energy, the confining potential energy and the inter-particle interaction 
energy, respectively. Here M and a denote the atomic mass and the s-wave scattering length, 
respectively, and i denotes the particle index. Let {|n)} be a complete set of exact eigenstates 
of H with eigenvalues {E n }, where n represents a complete set of quantum numbers that 
uniquely specify the many-body state of our system. For brevity of notation, we assume n 
to be nonnegative integers such that Eq < E\ < E2 < • • We are interested in low-lying 
excitations of many-body states which are excited by a general excitation operator F. Let 
|0) be the ground state and |1) be the lowest excited state with the excitation energy given 
by Tiuio = E\ — Eq. Because u w is not greater than huj n o = E n — E for n > 2, we obtain 
the following inequality 

2 < 2 \(MF\0)\ 2 + j: n ^\(n\F\0)\ 2 (^) 3 

\(l\F\0)\i + Z n ^\(n\F\0)\^ ' 1 j 

which shows that an upper bound hu) upper of the lowest excitation energy is given by 

hcu^ = y/m 3 / mi , (2) 

where m p = J2 n I (®\F\ n ) \ 2 (^ lu} no) p is the p-th moment of the excitation energyJll'S Similar 
methods which evaluate the excitation energy have been used extensively in the field of 
nuclear physics.0 The advantage of this formula is that m\ and m 3 can be expressed as 
expectation values of commutators between F and H with respect to the ground state |0) 
as mi = |(0|[Ft, [H,F}]\0) and m 3 = |<0|[[i^+, H], [H, [H,F]]]\Q), and we can therefore find 
hu uppeT without the necessity of finding excited stateslll To obtain u nppcr which is very close 
to co> 10 , we need to find a correct excitation operator F and an accurate ground-state wave 
function. 

We first consider the case of an axially symmetric trap (oj x = u y = ou±), which is 
relevant to recent experiments.iHi Fetter proposed a trial wave function for the condensate 
with repulsive interaction asil 



r 2 , Z 2 \ (l+r?)/2 



2 "2 



if jT^ + jTr2 < 1 and zero otherwise, where r± = ^a; 2 + I/ 2 , Co is the normalization constant, 

J__j_ 2 2 



= ^h/MtUj(j =_L, 2) is the oscillator length, and i?j and 77 are the variational parameters 
which are determined so as to minimize the total energy E = (T) + (U) + (H int ). Fetter 
showed that the trial wave function smoothly interpolates between the noninteracting 
limit (\I/(r) cx exp(— r 2 L /2d 2 L — z 2 L /2d 2 ,) for 77 ~ y/R±, \[R~ Z — » 00) and the strongly interacting 
limit (^(r) ~ c (l — r 2 jd\R\ — z 2 /d^R 2 ) 1 ^ 2 for 77 — > 0). However, how well it describes an 
intermediate regime had not been investigated until now. We therefore examine the accuracy 
of Fetter's variational wave function for various atomic numbers N . Table I compares 
the expectation values of the total, kinetic, potential, and interaction energies obtained by 
Fetter's wave function with those obtained numerically according to the method reported 
in ref. 28. The agreement is excellent and we thus conclude that eq. (3) serves our purpose 
very well. 

Experimentally, low-lying collective modes are excited by superimposing small ac cur- 
rents with appropriate phase relationship through magnetic coils. This leads to small mod- 
ulations of the frequencies of the confining potential. We consider \m\ = 2 and m = modes 
which are relevant to recent experiment sil~i The collective mode with magnetic quantum 
number \m\ — 2 describes the situation in which the condensate expands in one direction 
and simultaneously contracts in the other, thus maintaining its volume. This mode can be 
excited by modulating two radial trap frequencies out of phase but by the same amount 
Sou <ti u j_. The resulting perturbation which defines the excitation operator is given by 
F H=2 = M/2E, {[(wl + Su) 2 - u\] x 2 + [{u ± - 5uf - u 2 ± ] y 2 } « Muj^uj^x 2 - y 2 ). 
We note that F is proportional to Sj r « 2 (^2,2 + ^2,-2)) where Y/ m is the spherical harmonic 
functionl! We can therefore take the excitation operator as F\ m \ =2 = J2i( x j ~ yf), where 
the numerical factors are dropped because they are canceled in forming the ratio (2). The 
commutator [H, F| m | =2 ] is calculated as 
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[H, F\ m]=2 ] = — YJtPli + Pin x i ~ Vi] 



M 

We obtain the first moment m\ as 



2ih 

^(XiPxi-ViPyi)- (4) 



mi = ^([Fl ml=2 ,[H,F lm]=2 ] 

" / U*t + yt)) = T^r(u ± ), (5) 



where is a radial (x or y) component of the trap potential energy and (■ ■ •) denotes the 
expectation value over the trial wave function (0). We calculate the third moment m 3 as 

<(t±) + (u ± )), (0 



m 3 = 2\U^H=2>^J' n,Yn,r\ m \ =2 \ 



M 2 

where Tj_ is a radial (x or y) component of the kinetic energy. Because the upper-bound 
frequency does not directly depend on the interaction energy (H int ), we obtain the upper- 
bound frequency as 



w upper (m = 2)= u±y J 2 (l + (T ± )/(U ± )). (7) 

In the absence of inter-particle interaction we find that (T±) = (U±), so that u upper (m = 
2) = 2uj_, while in the TF limit we have (T±) = 0, so that eq. (|7|) reduces to the exact 
result cu uppcr (m = 2) = y/%} ± $ 

The m = mode describes an excitation in which the condensate alternately expands 
and contracts in the radial direction. This part of the excitation is described by F m=Q = 
M/2J2i{[( UJ ± + Sto) 2 — k>iJ r i.i} w Mcuj_5ujJ2i r ±i^ Because of the repulsive interaction, 
however, the condensate should also undergo oscillations in the axial direction which must 
be out of phase with the radial motion.! Hence we consider the excitation operator of the 
form F m=0 = Mu±8uJ2i( r ±i ~ az Y)-> where a is another variational parameter. We note 
that F is a linear combination of two modes J2i r i (n — 1,1 — m — 0) and J2i r iY2,o(9i, 4>i) 
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[n = 0, 1 = 2, m = 0). The excitation operator for the m = mode can therefore be taken 
as F m=0 = J2i( x i + yf ~ az i) f° r simplicity as in the \m\ = 2 case. The first moment mi is 
calculated as, 

m 1 = — h^L + a 2 ^), 8 

where U z is the axial component of the trap potential energy. The third moment is 
calculated as 

ms = ^ ([2(T± + 2U±) + q2{Tz + ^ + (1 " f )2Eint]) (9) 

where T 2 is the axial component of the kinetic energy. We thus obtain the upper-bound 
frequency as 



,u PP er, _ n . _ \ MT ± ) + (U ± )) + a\(T z ) + (U z )) + (1 - a/2f{H mt ) } \ 

[m-U,a)-[2 2(U±)/Lji + a?(U M )/u2 J' 1 Uj 

By minimizing u; upper (m = 0, a) with respect to a, we find u; upper (m = 0) = 2uj± in the 
noninteracting limit and cj uppcr (m = 0) = cj ± (2 + §A 2 - |V9A 4 - 16A 2 + 16)3 (A = uj z /uj±) 
in the strong interacting limit. This latter result is identical to that is obtained in ref.15 
by another method. The agreement shows that the excitation operator F m=0 that we chose 
is indeed correct. We also note that if uj npper is maximized with respect to a, we obtain 
uj uppcv (m = 0) = w_l(2 + |A 2 + |V9A 4 - 16A 2 + 16)1, which also coincides with the result 
reported in ref. 15 with higher frequency. Because states excited by F m= o are restricted to 
states which are constructed by linear combinations of the n = 1,1 — m — mode and the 
n = 0, 1 = 2, m = mode, there should be two values of a that make o; upper extremal and 
the corresponding states should describe the two lowest-energy excited states as obtained 
above. 

Our method can also be applied to the dipole (1 = 1) modest, which correspond to 
the center-of-mass motion of the condensate, and should therefore not be affected by the 
inter-particle interaction. This is known as the generalized Kohn's theorem. We consider 
the excitation operator F = M/2J2i UJ ±_[( x i + fi) 2 — X T\ ocJ^iXi'^- J2i r i(Yi,i( r i) — ^i,-i( r i)) or 



F = M/2£^i[(yi+<*) 2 -yf] ex ZiniYx^+Y^ri)) for \m\ = 1, and F = M/2£^[(^+ 
5) 2 — zf] oc £, z i Z)i r j^io( r j) f° r m = 0. Substituting F into Eq. (^]), we easily obtain 
w upper(j = 1; | m | = ^ = 0J± and yupper^ = l,m = 0) = u g . We thus find that the dipole- 
mode frequencies obtained by our method coincide with trap frequencies, being independent 
of the strength of interaction. 

Figure 1 compares our analytical results (solid curves) with the experimental data (dots) 
taken from ref. 5 for 87 Rb atoms, where we use the same parameters as in the experiment; 
a = 109a (a is the Bohr radius) and u z /y/8 = uo±_ = 2ir x 132Hz. Our results for both 
m = and \m\ = 2 modes are in excellent agreement with those of the experiment and 
with those obtained with the Bogoliubov approximatiorJl!. We have also calculated the 
upper-bound frequencies using numerically calculated values of (T) and (U) for 4500 atoms 
and obtained k; upper (|m| = 2) = 1A54uj± and u npper (m = 0) = 1.871a;j_. These results are 
in excellent agreement with those obtained using Fetter's variational wave function, that is, 
w upper(| m | = 2) = 1.450a; ± and w upper (m = 0) = 1.870w_i_. 

We briefly describe the results of our study for the case of a spherically symmetric trap 
(u>± = u z = uq), where we can compare our results with those reported in ref. 15 and with 
numerical ones based on the Bogoliubov approximation^. For the quadrupole mode (1 = 2), 
we consider the case of m = 2 without loss of generality, where F = M5uiJ2i( x2 ~ Ui)> an d 
find that u uppcT (l = 0) = 0JqJ2 + (T)/(U). For the monopole mode (n = 1, 1 = 0), where 
F = M5ujJ2irf, we find that uj uppci (l = 0) = oj ^5 - (T)/{U), where we have used the 
virial identity 2(T) —2(U) + 3(H int ) = 0. These results are the same as those reported in ref. 
15, although analytic evaluation of (T) and (U) in an intermediate strength of interaction 
is not so far available from the TF limit. 

In Fig. 2, we present our results for both modes with the same parameters (uq = 
2tt x 200Hz, a = llOao) as those reported in ref. 17 using the Bogoliubov approximation. Our 
results are almost indistinguishable from those obtained numerically using the Bogoliubov 
approximation. 

In conclusion, we have presented an analytical method to evaluate almost exactly the 

6 



collective-mode frequencies of a trapped Bose gas. The results are in excellent agreement 
with those of the Bogoliubov approximation!!! and those of the experiments!! at zero tem- 
perature. 
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TABLES 

TABLE I. Comparison of total energy (E to t), kinetic energy (T), confining-potential energy 
(U) and interaction energy (iiint) obtained using Fetter's variational wave function with those 
obtained numerically according to ref. 28 (in parentheses) for an axially-symmetric trap. We take 
y/8uj± = u) z = 2-7T x 220 Hz and a = lOOao- Energies are shown in units of huj± and No is the atom 

number. 

N (E tot ) (T) (U) (H int ) 

1000 3.86 (3.84) 0.75 (0.76) 2.17 (2.15) 0.95 (0.93) 

10000 7.83 (7.76) 0.46 (0.45) 4.61 (4.57) 2.76 (2.74) 

20000 10.06 (9.98) 0.40 (0.38) 5.95 (5.91) 3.70 (3.68) 
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FIGURES 

FIG. 1. Collective- mode frequencies of the \m\ = 2 and the m = modes in an axially 
symmetric trap, where we assume w 2 /V8 = u>± = 2ir x 132Hz and a = 109ao- The solid curves 
show our results, the dots show the experimental data reported in ref. 5, and the dashed lines 
show the results in ref. 15. 

FIG. 2. Collective-mode frequencies for the quadrupole and the monopole modes in a spheri- 
cally symmetric trap. The dashed lines show the results reported in ref. 17. 
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